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1. INTRODUCTION 

The strong coupling as is a fundamental parameter of the Standard Model. 
In comparison to parameters like aom, Mz and sin^ 9w it is relatively poorly 
known. However the precision of as measurements has improved dramatically 
in recent years. More than twenty different types of process, from lattice QCD 
studies to the highest energy colliders, can be used to measure as accurately. 
The most precise determinations now quote uncertainties in as{M'^) of less 
than 5%. There is also a remarkable consistency between the various measure- 
ments. 

A comprehensive review of as measurements, including detailed descriptions 
of the underlying physics for the most important processes, can be found in 
Rcf. One year later, several of the measurements quoted in Ref. [Q have 
been updated, resulting in a slight shift in the overall 'world average' value. The 
purpose of the present review is to update the discussion on as measurements 
given in Ref. Q, focusing on the new values reported in the last year. For 
more theoretical details, descriptions of other measurements and a full set of 
references, the reader is referred to the original review in Ref. 

The current situation is summarised in Fig. |l|, which updates Table 12.1 of 
Rcf. ||l|. Before discussing the new measurements in detail, we begin with some 



(* ) Based on a talk presented at the 'New Non-Perturbative Methods and Quanti- 
sation on the Light Cone' conference, Les Houches, France, March 1997 
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Fig. 1. 
Ref. [h 



Measurements of Qg(Mf ), in the MS renormalisation scheme, updated from 



technical preliminaries. In perturbative QCD the dependence of the strong 
coupling on the renormalisation scale is determined by the /3-function: 



Q 



2das{Q^) 



aQ2 
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(3{as) - -bal.{l + b'as + b"al + ...) , (1) 

where b = (33 — 2nf) / (127r) etc. The coefficients in the perturbative expansion 
depend, in general, on the renormahsation scheme (RS), although for massless 
quarks the first two coefficients, b and b' , are RS independent. In essentially 
all phenomenological applications the MS RS is used; see Ref. for further 
discussion and explicit expressions for the known /^-function coefficients. 

At leading order, i.e. retaining only the coefficient b, Eq. (|l|) can be solved 
for as to give 

"^^"^^ l + as{Ql)bHQyQl) 

or 

= bhWWr 

These two expressions are entirely equivalent - they differ only in the choice of 
boundary condition for the differential equation, as{Qo) in the ffist case and 
the dimensionful parameter A in the second. In fact nowadays A is disfavoured 
as the fundamental parameter of QCD, since its definition is not unique beyond 
leading order (see below), and its value depends on the number of 'active' quark 
flavours. Instead, it has become conventional to use the value of as in the MS 
scheme at = as the fundamental parameter. The advantage of using 
Mz as the reference scale is that it is (a) very precisely measured Q , (b) safely 
in the perturbative regime, i.e. as{M'^) 1, and (c) far from quark thresholds, 
i.e. rrib < Mz < mt- 

The parameter A is, however, sometimes still used as a book-keeping device. 
At next-to-leading order there are two definitions of A which are widely used 
in the literature: 

definition 1 : bln^ = + b' In ^ ^'"f^^i,, ) , (4) 

A2 as{Q^) \l + b'as{Q^)J 



definition 2 : as{Q^) 



1 



51n(g2/A2 



b' lnln(QVA2^ 



b In(g2/A2 



(5) 



The first of these solves Eq. (|l|) exactly when b" and higher coefficients are 
neglected, while the second (the 'PDG' definition j2j) provides an explicit ex- 
pression for as{Q^) in terms of Q"^/ A"^ and is a solution of Eq. (|^) up to terms 
of order 1/ ln^(Q2/A2). Note that these two A parameters are different for the 
same value of as{M'^), the difference being about one quarter the size of the 
current measurement uncertainty: 

^(5) _ ^(5) _ _ i^^^pA(^). (6) 

In this review we will be mainly concerned with measurements from e^e^ 
colliders (in practice LEP and SLC) and from deep inelastic scattering. Both 
processes offer several essentially independent measurements, summarised in 
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Table I. — Summary of the most important processes for as determinations in e e 
collisions and in deep inelastic lepton-hadron scattering. 





quantity 


perturbation series 




Ree, Rz, Rt 

event shapes, /a, ... 


R = Ro[l + as/7T + ...] 
1/ada/dX = Aas + -Ba| + . . . 
dD''/d In Q2 = asD'' P + . . . 


£N DIS 


cr(2 + 1 jet) 


dF,/d\nQ^ ^asF,(E)P+ ... 
J dxF,{x, Q'^) = A + Bas + ... 
a = Aas + Bag + . . . 



Table ffl. Note that aU of these use the qqg vertex to measure as, with the 
high scale provided by an electroweak gauge boson, for example a highly 
virtual 7* in DIS or an on-shell boson at LEPl and SLC. There are two 
main theoretical issues which affect these determinations. The first is the ef- 
fect of unknown higher-order (next-to-next-to-leading order (NNLO) in most 
cases) perturbative corrections, which leads to a non-negligible renormalisation 
scheme dependence uncertainty in the extracted as values. This is particularly 
true for the event shape measurements at e'^e^ colliders. The exceptions here 
are the total e+e^ hadronic cross section (equivalently, the hadronic decay 
width) and the DIS sum rules, which are known to NNLO. The second issue 
concerns the residual impact of ©(l/Q") power corrections. For some processes 
it can be shown that the leading corrections are 0{1/Q) (for example 0{1/Mz) 
for the corrections to event shapes at LEPl and SLC) which can easily be com- 
parable in magnitude to the NLO perturbative contributions. In deep inelastic 
scattering, the higher-twist power corrections are 0(1/Q^(1 — x)) and must be 
included in scaling violation fits especially at large x. Such power corrections 
(and their uncertainties) must be taken into account in as determinations, 
either using phcnomenological parametrisations or theoretical models. 

Before discussing the new high-energy collider measurements of as it is im- 
portant to mention also determinations from lattice QCD, which have very 
small uncertainties. One of the simplest ways to define as on the lattice is to 
use the average value of the 1x1 Wilson loop (plaquette) operator: 

An /3.4\ 

lnWi,i = — ap — [1- (1.19-H0.07n/)ap], (7) 
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where a is the lattice spacing. A variety of choices is available for determining 
a, i.e. measuring the scale at which ap has the value measured in Quarko- 
nium level splittings, for example T{S ~ P) and T(15' — 25"), are particularly 
suitable. Subsequently the plaquette ap can be converted to the standard MS 
as for comparison with other determinations: 



a 



^MS,„/)(Q2)^^(rv)(g5/3g2) 



1 



-a 



(«/) 



C2{nf) (c 



(8) 



At present the two-loop coefficient is known only for = - the shift in as 
between using 6*2(^1/ = 0) and C2 = can be used to define a 'conversion' error. 
Several new lattice as values have been obtained recently, see for example 
Ref. [Q , and are included in Fig. |l|. As an example of the high precision of these 
measurements, we quote the value obtained by the NRQCD collaboration |^ 
using the T(5 — P) splitting: 

af^(M|) = 0.1175 ± 0.0011(stat. + sys.) ± 0.0013(m^y" ) ± 0.0019(conv.), (9) 

where the first error is due to the lattice statistics and systematics, the second 
is from the extrapolation in the dynamical quark mass, and the third is the 
conversion error mentioned above. 

In the following sections we will discuss new as measurements from LEP/SLC 
and from deep inelastic scattering. Section ^ presents a new value for the as 
world average. 



2. as from LEP and SLD 



In principle the most reliable determination of as at the LEP and SLD e^e 



colliders comes from the Z° 
ratio Rz, 



hadronic width. In particular we have, for the 



Rz 



r(Z" ^ hads.) 



Ro 



as 

TT 



m 



(10) 



with Ro = SEqivl + al)/{v. 



aj). The perturbative coefficients are known 
up to third order, see Ref. [y for explicit expressions and references, and as a 
result the prediction is very stable with respect to variations in the renormali- 
sation scale. In practice, since Rq depends on the weak mixing angle and other 
electroweak parameters, it is more appropriate to perform a global fit to all rel- 
evant electroweak quantities, for example the (LEP and SLD) partial widths 
and decay asymmetries, pp collider measurements of Mw and rrit, etc. Such 
analyses are performed regularly by the LEP Electroweak Working Group, and 
the results of a recent (1996) fit 1^ are summarised in Table ||. An additional 
theory error from unknown higher-order corrections of das{M^) = ±0.002 has 
been estimated, see Ref. B and references therein. The resulting as value. 



as(Af|) = 0.120 ± 0.003(fit) ± 0.002(theory), 



(11) 
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Table II. — Values for the Standard Model parameters obtained from a global fit to 
LEP, SLD, pp and uN data, from Ref. g. 



parameter 


fit value 


mt [GeV] 
Mh [GeV] 
«s(M|) 


172 ±6 
0.120 ±0.003 


sm ej 
1 - m2,/m| 

Mw [GeV] 


0.23167 ± 0.00023 
0.2235 ±0.0006 
80.352 ±0.033 



is displayed in Fig. 

The other high-precision determination of as at LEP and SLC comes from 
event shapes, quantities which measure the relative contribution of the 0{as) 
e^e^ — !■ qqg process to the total hadronic cross section, see Table 0. A typical 
example is the thrust distribution: 

i ^ = asA^iT) + alA,{T) + ... + (^-i-) . (12) 

Such quantities are known in perturbation theory to ©(a^), and the theoretical 
predictions in the T I region can be improved by resumming the leading 
logarithmic An ~ ln^^"^^-'(l — T)/(l — T) contributions to all orders, as dis- 
cussed in Ref. (l). Another important recent theoretical development has been 
an improved understanding of the leading 0{1/E) power corrections Q, which 
at LEP can be as numerically important as the next-to-leading perturbative 
corrections. 

Event shapes have yielded as measurements over a wide range of e+e^ 
collision energies, the most recent measurements being at the LEP2 energies 
^/s = 161 and 172 GeV. Although the statistical precision of these measure- 
ments cannot match that obtained at the Z*^ pole, the results are consistent 
with the evolution of as predicted by Eq. (|l|). For example, Fig. ^ shows 
the as values determined by the L3 collaboration ||] from event shape mea- 
surements at LEPl and LEP2 energies. The solid line is the evolution predicted 
by perturbative QCD. Figure |l| contains a new 'LEPl. 5' average value for as 
obtained from event shapes at y/s = 133 GeV, taken from the 1996 review by 
Schmelhng §: 

as(Q2 ^ (133 GeV)2) = 0.114 ±0.007 ^ as(A/|) = 0.121 ± 0.008. (13) 

Another updated value in Fig. |l| is that obtained from the scaling violations of 
the fragmentation function measured in e+e^ hX over a range of collision 
energies, the analogue of the scaling violations of structure functions in DIS. 
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The new value (an ALEPH/DELPHI average taken from Ref. corresponds 
to 

as(M|) = 0.124 ±0.010. (14) 

Finally, the CLEO collaboration have published jl^ a new value for as 
obtained from the relative decay rate of the T(15) into a single hard photon: 



p799 



a 



1- (2.6-2.11n(mg^2)) 



Psff 5 asifJ.'^) 
The new value, 

as{M^(is)) = 0.163 ± 0.002(stat.) ±0.014(sys.) 

^ as{M^) = 0.110 ±0.001(stat.)±0.007(sys.), 
is included in Fig. [l]. 



(15) 



(16) 




Fig. 2. — Measurements of as from event shapes at LEPl and LEP2 from the L3 
collaboration M. The errors correspond to experimental uncertainties. 



3. as from Deep Inelastic Scattering 

The traditional method of measuring as in deep inelastic scattering is from the 
strength of the structure function scaling violations predicted by the DGLAP 
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equations: 

Q'-^ = ^^^{P^'^^q' + P^^^^g), (17) 

where g^"^ and g'^ are respectively non-singlet and singlet combinations of 
quark distribution functions. The fixed target and HERA structure function 
data, spanning a large range in x and Q^, are all consistent with NLO DGLAP 
evolution, and yield 015 values which are in broad agreement. As an example, 
Fig.|| ^] shows the values for various DIS data sets as a function of the 
as{M^) value in the evolution equations. With one exception, all data sets 
exhibit a minimum in the as — 0.11 — 0.13 range. In fact the 'best fit' value 
for these data sets is as{M'^) = 0.118, exactly the world average value (see 
Section || below). 

It is difficult to extract a proper error on from such global fit analyses. 
This requires a rigorous treatment of systematic errors and inclusion of higher- 
twist contributions in the fit. Several groups have performed such analyses. 
For example, the Milsztajn-Virchaux analysis of the SLAC/BCDMS (eiV, iiN) 
data |l^ yields 

as(M|) = 0.113 ±0.005, (18) 

where the error includes statistical, systematic and scale dependence uncer- 
tainties. Recently the CCFR collaboration have reported ||l^ a new value of 
as from their Fj'^, xF^^ high-precision data (see Fig. |^): 

as(M|) = 0.119 ± 0.002(exp.) ± O.OOl(HT) ± 0.004(scale). (19) 

The second error is from an estimate of the higher-twist contribution using the 
model of Ref. iQ, and the third is the scale dependence uncertainty imple- 
mented as in Ref. ||l2| . Note that the value in (|9|) is somewhat larger than the 
earher (1993) CCFR value of as = 0.111 ± 0.004. The change is due to new 
energy calibrations of the detector |l^ . 

Deep inelastic scattering structure functions satisfy a variety of sum rules, 
corresponding to the conservation of various nucleon quantum numbers. In 
general the parton model values of the sums have ©(as) corrections, which can 
be used to extract as from measurements of structure function integrals at fixed 
Q^. Two sums rules which have been used to obtain precision measurements 
are the Gross-Llewellyn Smith and Bjorken sum rules (see Ref. ^ for more 
discussion and references): 



GLS: / dx{F^P + F^P) = 6 1 + -^ + ... + Aht, (20) 







as 
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Fig. 3. — values for various DIS data sets obtained in a global fit to these and 
other hard scattering data jl^] . 



BjS: 



&9V 



as 



iHT, 



(21) 



where Aht represents 0(1/(3^) higher-twist contributions. A new analysis 
of polarised structure function measurements has produced an update of the 
as value from the Bjorken sum rule. In Ref. Padc Summation is used 
to reduce the theoretical error from the choice of renormalisation scheme in 



the calculation of the perturbation series on the right-hand side of (21). The 
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Fig. 4. — Measurements of the structure function F2 from the CCFR collaboration 
together with a NLO QCD fit, from Ref. 



resulting theoretical error in as{M'^) is estimated at ±0.002: 

as(M|)=0.117t^; 
This new value is included in Fig. 



as(M|) = 0.117t^ ™^(exp.) ± 0.002 (theory). (22) 
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Finally, as can be obtained from jet fractions and event shapes in DIS, see 
Table |l|. For example, NLO theoretical predictions are currently being used at 
HERA to extract as from the relative rate of '2+1' jet production at high Q"^, 
the analogue of /a in e+e^ annihilation. No new results have been published 
since the review in Ref. 

4. Summary 

The average valu^Ql of the measurements presented in Fig. |^ is 

WORLD AVERAGE: as(M|) = 0.118 ± 0.004. (23) 

Following Ref. [Q, the error here is defined as 'the uncertainty equal to that of 
a typical measurement by a reliable method'. In view of the recent improve- 
ments in the lattice, Z° hadronic width, and DIS {i^N) determinations, it seems 
appropriate to decrease the uncertainty of ±0.005 in Ref. to ±0.004. The 
central value in (^3|) has increased by +0.002 from that given in Ref. §. This IS 
due primarily to (a) increases of +0.003 and +0.004 in the central values of the 
two lattice determinations, and (b) an increase in the CCFR i/N DIS scaling 
violation central value of +0.008. In view of the remarkable consistency of all 
the measurements, and in particular of those with the smallest uncertainties, it 
seems unlikely that future 'world average' values of as will deviate significantly, 
if at all, from the current value given in (^). 
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